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Abstract 



New results are presented for the degeneracy condition of elastic waves in 
anisotropic materials. The condition for the existence of acoustic axes involves a 
traceless symmetric third order tensor that must vanish identically. It is shown 
that all previous representations of the degeneracy condition follow from this 
C^ , acoustic axis tensor. The conditions for existence of acoustic axes in elastic crys- 

tals of orthorhombic, tetragonal, hexagonal and cubic (RTHC) symmetry are rein- 
terpreted using the geometrical methods developed here. Application to weakly 
^ . anisotropic solids is discussed, and it is shown that the satisfaction of the acoustic 

O I axes conditions to first order in anisotropy does not in general coincide with true 

acoustic axes. 



cn ■ 1 Introduction 
in 

^ I An acoustic axis is defined as a direction in which two sheets of the slowness surface 

O ' intersect. Propagation directions for which this occurs permit a pair of degenerate po- 



larization vectors associated with the common slowness, and hence circularly polarized 
waves may propagate. The importance of acoustic axes in crystals has been long appre- 
ciated, and the conditions required for their existence has been known since the early 
1960's. In a paper on this subject in 1962, Khatkevich |1] derived conditions for the 

'^ ' existence of acoustic axes and discussed their application to the different crystal symme- 

try classes. Acoustic axes are "special directions" in the sense of Fedorov |2] , similar to 
O ■ directions in which transversely polarized waves propagate. In fact, an acoustic axis is 

a direction for transverse waves, but not vice versa in general. There has been renewed 

'k>( I interest in acoustic axes in the past decade, initiated in part by Fedorov and Fedorov 

Vh I [S] who derived new and simplified representations of the acoustic tensor. Based on this 

approach, Boulanger and Hayes jH Ej provided a complete classification of the acous- 
tic axes in elastic crystals of orthorhombic, tetragonal, hexagonal and cubic (RTHC) 
symmetry. In an interesting but unrelated paper, Mozhaev et al. p] gave an explicit 
condition for crystals of trigonal symmetry, thereby completing the classification for all 
crystal symmetry classes higher than monoclinic. 

The purpose here is not to attempt to extend this classification, but to focus on the 
general conditions and how they are formulated. In particular we provide a definitive 
set of conditions which must be satisfied for acoustic axes to exist in an elastic material. 
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A note on notation: lower case boldfaced latin symbols denote vectors, and vectors 
are normally of unit length, such as the unit direction vector n, and the orthonormal 
triad {61,62,63} which define the fixed directions. 

There are several ways of representing a a second order symmetric tensor A. We 
begin with the perhaps less familiar Hamilton's cyclic form [7j, also known as Hamilton's 
decomposition jS] and the biaxial form j2]: 

A = AI + — (ai®a2 + a2®ai). (1) 

The numbers A, A', and the biaxial unit vectors ai,a2 are uniquely defined by A. The 
proof of eq. (Hj) may be seen by direct construction, see Appendix A. Fedorov calls 
the directions ai and a2 the axes of the tensor (not to be confused with the principal 
axes). The standard diagonal form for A, also known as the Gibbs decomposition [Hj, 
employs the principal axes, an orthonormal triad {a, b,c}, 

A = Aia(g)a + A2b(g)b + AscOc, (2) 

with the eigenvalues ordered Ai > A2 > A3. Equations ^ and (j21) can be reconciled 
using the identity a®a + b(8)b + c®c = I, implying A = A2, A' = Ai — A3, and 

'^J-iT^r-'^+i-m^f'c, , = 1,2. (3) 

Ai — A3 Ai — A3 

Any positive definite symmetric tensor A defines an ellipsoidal surface {x : x • Ax = 
1}. A symmetric second order tensor is uniaxial if two of its eigenvalues coincide and the 
third corresponds to the preferred axial direction. The ellipsoid reduces to a spheroid, 
which may be oblate or prolate, depending as Ai = A2 > A3 (oblate) or Ai > A2 = A3 
(prolate). In either event, A reduces to 

A = A2I ± A'd (g) d, uniaxial, (4) 

where, according to eq. (jT}, the ± corresponds to ai = ±a2, and the implications for 
the eigenvalues Ai, A2, and A3 and the definition of d follow from eq. (jH)). Note that the 
sign of the d ® d term in eq. (0} directly indicates whether the uniaxial tensor is prolate 
(+1) or oblate (—1). In the limit as the three eigenvalues coincide eqs. ((H), Q and (JH) 
all reduce to A = A2I as A' -^ 0, i.e. A is spherical. An oblate (prolate) acoustical tensor 
corresponds to one where the two larger (lesser) phase speeds coincide for a given choice 
of the phase direction. It is worth mentioning that the geometric descriptions (oblate, 
etc.) do not refer to the slowness surface but to the form of the acoustical tensor for the 
given phase direction. 

Our focus will be on discussing the conditions under which the acoustic tensor of 
elasticity, Q defined below, is uniaxial. Before doing so, it is useful to note that it is 
relatively simple to define conditions for sphericity. Thus, the identity 

(Ai - A2)' + (A2 - A3)' + (A3 - Xif = 2 /' - 6 // , (5) 
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where / and II are invariants of A, J = tr A, 211 = (trA)^ — trA , implies that 
all three eigenvalues coincide iff the condition P = 311 is met. Quantities such as 
(Ai - Aa)^ + (A2 - As)^ + (A3 - Ai)3 and T = (Ai - A2)(A2 - A3) (A3 - Ai) vanish iff A is 
uniaxial. However, it does not seem possible to represent these in terms of the invariants 
of a second order tensor: /, // and the third scalar invariant III = detA. We will find 
instead that a third order tensor is required. 

2 Acoustic axes 

The acoustical tensor of elasticity, Q(n), is symmetric positive definite and defined as 
Qjfc(n) = CijkiUjni, where C with elements Cijki is the elastic stiffness. The eigenvalues 
of Q are Xj = pv^, where p is the mass density and Vj{n), j = 1, 2, 3 are the phase speeds 
for phase direction n. According to the general representation of eq. Q, at an acoustic 
axis, Q must be of the form 



(6) 



where we have chosen for purposes of clarity to write the six elements of Q as a column 
"vector" . The form of Q in eq. (jHl) is motivated by the exhaustive study of Boulanger 
and Hayes jH Ej for RTHC crystals. Boulanger and Hayes, using general expressions 
derived by Fedorov and Fedorov |3j, wrote Q as in eq. ^ and obtained explicit results 
for the number and orientation of acoustic axes in RTHC crystals. We do not restrict the 
present analysis to any particular elastic symmetry, although we will return to the RTHC 
crystals later in Sections and IHl Focusing on the final three elements in Q, implies 
that q2q3 = Q23, etc. For the moment let us assume that all three of Q23, Q31, Q12 are 
non-zero. Define the vector q(n) = giei + ^2^2 + ^'363 by 



' Qii ■ 




' 1 ~ 




r 9? 1 


Q22 




1 




qI 


^ ^33 
Q23 


> = X2 < 


1 




> ± < 


^3 . 
^293 


Q3I 









^391 


[ Q12 J 









. 9192 , 



Q2 
93 



Q31Q 



12U/2 



^23 ' 

Ql2Q23^|l/2 

Q 



signQ23, 



31 



Q23Q31n1/2„. 
Q12 ^ 



signQsi, 
sign(5i2. 



(7) 



where s = sign(523Q3iQi2- The acoustic tensor then can be written. 



Q(n) = D + sq®q, 
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where D is diagonal with Dj = 

Q31Q12 



Qii ~ sqf (no sum), that is, 



Di = Q 



11 



Q 



D2 = Q 



Q12Q 



23 



22 



23 
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D, 



Q33- 



^23^ 



31 



31 



Q 



(9) 



12 



The expression of the acoustical tensor in the form (jHJ, while entirely general only as 
long as Q23Q31Q12 7^ 0, is very useful for examining at this stage some, but not all, of 
the conditions under which acoustic axes occur. 

In particular, Q is uniaxial when D is spherical, that is: Di = D2 = -D3. If this 
occurs, then the uniaxial acoustical tensor is prolate or oblate depending as s = +1 or 
s = —1, respectively. This implies that the two outer (inner) sheets of the slowness 
surface coalesce if s = +1 (s = —1). Furthermore, the vector q is in the direction of 
the third polarization, perpendicular to the plane of polarization of circularly polarized 
waves. Equating the Di implies the conditions, 



D, 



D-, 



Di 



D3 


=^ 


Q22-Q33 -- 


= Q23 


Di 


=^ 


Q33-Q11 -- 


= Q31 


D2 


=^ 


Qii — Q22 = 


= Q12 



.<5l2 Q31- 



'Q31 

'Q23 
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Q12' 
Q12, 
Q23' 
^23. 



(10) 



Q23 Q31' 



In order to account for the possibility that the denominator vanishes, these equations 
are normally expressed [Q 

0, = 0, J = 1,2,3, (11) 

where (pj are found by multiplying eqs. (fTUj) by Q12Q23Q31, 



Q31Q12 {Q33 
Q12Q23 [Qii 



Q22) 
Q33) 



Q23 {Q12 ~ <53l)' 
^31 {Q23 ~ Q12)' 



(12) 



Q23Q31 {Q22 — Qii) + Q12 {Q31 ~ Q 



23 J 



Equations (|lip were first derived by Khatkevich in 1962 ^ by a different method. 
Khatkevich argued that the orientation of acoustic axes in general follows from the 
condition that the cofactor matrix of M = Q — AI must vanish. Only three of the 
elements of the cofactor matrix are independent because of the simultaneous requirement 
that det M = 0. Setting the off-diagonal elements of the cofactor matrix to zero implies 



0, 



(13) 



Qn — A Q12 _ Q12 Q22 — A _ Q12 Q23 
Q13 Q23 Q13 Q23 Qi3 Q33 ~ ^ 

and eliminating A from these recovers Khatkevich's conditions (fTT|) . There has been 
some confusion about the starting point for Khatkevich's derivation, i.e., the fact that 
cofM = at an acoustic axis. Mozhaev et al. jH] note that this assumption is not 
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proved explicitly by Khatkevich or in the references that he offered, and they suggest 
that the arguments given by Fedorov for this condition are "rather tangled". However, 
at an acoustic axis Q has the form of (jHI with D = AI. Therefore M = s q q is of 
rank one, and the cofactor matrix is identically zero: cof Mij = eikiejmnQkQiQmQn = 0. 

It has been argued P that only two of the three conditions (fTTj) are independent, 
since, for example, the third equality Di = D2 follows from the previous two in eqs. 
(|12|). However, the conditions (jllj) are not always sufficient, as for instance, on a plane 
of symmetry. Consider n lying in a material plane of symmetry: to be specific, suppose 
63 is normal to a plane of material symmetry and ns = 0. Then Q23 and Qsi both vanish, 
and the three conditions ()1H) are identically satisfied. An additional equation is clearly 
necessary. One may be constructed by taking the product of the first two identities in 
eq. (Uni): 

(Qii - Q33){Q22 - Qss) = QI2 (1 - S)(l - S)- (14) 

V12 'V12 

If the elastic body possesses a plane of symmetry with normal e, and if n lies in the 

plane of symmetry, then Q(n) must display the same symmetry, viz. m • Qe = for all 

m ± e. With e = 63, and ^3 = 0, eq. (fT^ becomes 

(Qii - Q33) (Q22 - Qsa) - QI2 = 0. (15) 

This condition was also derived by Khatkevich pp, using a different argument. We show 
later that it is a consequence of another condition (see eq. (jS^ ). and will return to it 
after a more consistent and thorough set of conditions have been discussed. The point 
at this stage is to emphasize the multiplicity of conditions, some for general directions, 
others for planes of symmetry. This apparent complexity raises the question of what 
are the minimum number of conditions, and how may they be expressed? One solution 
was offered by Al'shits and Lothe [H] who derived a set of conditions in terms of a seven 
component vector which must vanish for acoustic axes to occur in a given direction n. 
An alternative solution to the question will be provided next, and the Al'shits and Lothe 
conditions will be considered later and shown to be overly restrictive. 

3 A coordinate invariant set of conditions for acous- 
tic axes 

Consider a symmetric tensor A with orthonormal eigenvectors a, b and c, and distinct 
eigenvalues Aj, 2 = 1,2, 3. Noting that 

A^ = Aia(g)a + A2b(g)b + A3C(g)c, (16) 

implies, for instance, that Am and A^m G span{a, b} for any vector m perpendicular 
to c. The triple product of three vectors, [p, q, r] = eijkPiqjTk, is zero if the vectors 



Acoustic axes in elasticity 6 

are coplanar. Hence, if m G span{a, b}, it follows that [m, Am, A m] is zero. Conse- 
quently, the vanishing of 

[m. Am, A^m] (17) 

is a necessary consequence of the fact that m is perpendicular to one of the eigenvectors 
of A. It remains to show that it is also a sufficient condition. 

In order to understand the function of m in eq. (|T7jl first consider the more general 
quantity [m, Ap, A q] = m • (Ap A A q) for arbitrary m, p and q. Direct calculation 
based on eq. (jlfjj) . and the fact that {a, b, c} form an orthonormal triad, gives 

Ap A A^q = (Asb • pc • q — A2C • pb • q) A2A3a^- 
(Al c • p a • q - A3 a • p c • q) AsAib + 
(A2a ■ pb ■ q — Ai b ■ pa ■ q) A1A2C. (18) 

The general form of (fTTjl follows from the identity (|TH|l as 

[m. Am, A'm] = (Ai - A2)(A2 - A3)(A3 - Ai) (m ■ a) (m ■ b) (m ■ c) . (19) 

By assumption, the eigenvalues are distinct, implying that the factor premultiplying 
(m • a) (m ■ b) (m ■ c) is non-zero. Therefore, if the triple product in eq. ()19|) vanishes 
the vector m must be perpendicular to at least one of a, b or c. This proves that the 
triple product (fTTj) vanishes iff m is orthogonal to one of the eigenvectors of A when the 
eigenvalues are distinct. 

Now consider the special case that A is uniaxial, say A2 = A3. Then every direction 
perpendicular to a is an eigenvector: this is circular polarization corresponding to the 
circle of eigenvectors {p : |p| = 1, p • a = 0}. Conversely, for a given arbitrary direction 
m, there is an eigenvector on this circle that is perpendicular to m, specifically p || mAa. 
Hence, [m. Am, A^m] vanishes for any and all directions m if A is uniaxial. The 
converse follows from the explicit form ()19|) . Thus, if the triple product vanishes for all 
m then 

(Ai-A2)(A2-A3)(A3-Ai) = 0, (20) 

which implies that at least two of the eigenvalues are equal. This proves the fundamental 

result: 

Uniaxial Condition: The symmetric tensor A is uniaxial if and only if the triple 

product [m. Am, A^m] vanishes for all m. 

This can be immediately applied to the acoustic tensor Q(n), yielding the following 
criterion for the existence of an acoustic axis: 

n is an acoustic axis <^=^ /(m, Q(ii)) = O5 V m, (21) 

where the function / is defined 

/(m) = /(m, Q(n)) = [m, Qm, Q^m]. (22) 
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Referring to eq. p9|l we see that this is in general an overly stringent condition. Specif- 
ically, it need not apply to all directions m but only to a sufficient number of directions 
chosen so that (m ■ a) (m ■ b) (m ■ c) 7^ 0. This requires that m not be limited to the 
set {a, b, c}, which is satisfied if, for instance, eq. (PT|) holds for four m, as discussed 
in Section m 

4 The acoustic axis tensor 

The function /(m, Q(n)) is third order in m and defines a third order tensor 0(n) via 
eqs. (Uni) and (ES): 

= (Ai - A2)(A2 - A3)(A3 - Ai) {abc} , (23) 

where {abc} is the totally symmetric part [lOj of a gD b C?) c, i.e. 

{abc},^, = I ia.b,c, + aM^ + aAc, + a^c, + a^C + a,h,c.). (24) 

In terms of the fixed coordinate system, 

/(m) = (pijk miirLjink, and = (pijk ^i ® e^- (g) e^ . (25) 

The tensor </> is totally symmetric, that is, (pijk is unchanged under any permutation 
of the suffices i, j, k and hence has 10 distinct elements ipijk- However, from its definition 
in eq. ((23), 

<Pikk = for z = 1, 2, 3. (26) 

A totally symmetric tensor with this traceless property is called harmonic ^lOj by virtue 
of the fact that V'^cpijkXiXjXk = 0, which follows from eq. (PB|) . The three relations (PB|1 
imply that (j) is characterized by seven independent elements. These may be expressed 
via an explicit decomposition of Dew, the space of symmetric traceless third order tensor. 
This is achieved by a partition Def into a direct sum of spaces of one or two dimensions 
each of which is irreducible under rotation about an axis, i.e. the group SO (2). The 
procedure, called "Cartan decomposition" when applied to harmonic polynomials, has 
been found useful in describing the symmetry classes of elasticity ^T] and photoelasticity 



Specifically, Dev is the direct sum of Dq, Di, D2 and D3, the spaces spanned, respec- 
tively, by Eq, and the pairs (Ei, E2), (E3, E4) and (E5, Eg), where 
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= (3/20)^ 


/M(5er 
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The seven spaces are grouped into the subspaces D„, n = 0,1,2,3, because of the 
property that elements of the associated spaces are rotated through angle n6 under a 
coordinate rotation of 6 about the axis 63. The reader is referred to Appendix B for a 
derivation of this result. Furthermore, the tensors E„ form an orthonormal basis of Biev, 
and hence the general representation of the acoustic axis tensor, or any harmonic third 
order tensor, is 

6 

(/) = ^ c„ E„ where c„ = E„ ■ . (28) 

n=0 

Direct calculation using eqs. PH|1 and (J2ZI) yields 

Co =(5/2)1/2 0333, Ci = (15/4)1/2 0133, 

C2 = (15/4)1/2 0233, C3 = (3/2)1/^ (0^3 - 0223), (29) 

C4 = 6^/2 0123, C5 = (0111 - 30i22)/2, 
C6 = (0222 - 30ii2)/2. 

Inner products of two tensors are cf) ■ (p' = (f)ijk(pijk = cqCq + Cic'i + . . . + C7C7. The strict 
condition for acoustical degeneracy, = is therefore equivalent to the requirement 
that the seven coefficients Cj, j = 0, 1, 2, . . . , 7 vanish ^. This also follows from the 
general identity 

6 

101' = E^- (30) 

n=0 

The elements (p may be expressed in terms of the elements of the acoustical tensor 
Q(n). Direct calculation using eq. (j^ gives 

0m = /(ej) = 0i (no sum), (31) 

where 0,, i = 1,2, 3, are defined in eq. (fT^ . Similarly, it can be shown that 

30nj = [{Qjj-Qii)(l)i-Qij(pj + Qik(pk]/Qij, i^j^k=^i (no sum), 

(32) 

60123 = (Q11Q22 - Ql2)iQ22 - Qu) + (^22^33 ' ^23)^33 " Q22) 
+ (Q33Qll-Q3l)(Qll-Q33). 

For instance, 0i23 is the only nonzero element of if Q is diagonal. Also, if Q23 and 
Qai are zero, as is the case when the material has a plane of symmetry with normal 63 
and n _L 63, then 

0aa3 = -(-1)" [{Qu - Q33)(Q22 - Q33) " QL1Q12/3, « = 1 Or 2. (33) 

Hence, the condition (J15p . which was derived in an ad hoc manner, can be viewed as a 
consequence of 0ii3 = or 0223 = 0. 
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The elements (pijk are related to the seven quantities Ri, i = 1, 2, . . . , 7 of Al'shits 
and Lothe fH^ as follows: 

Rl = —0333, R2 = 0222, -R3 = —0111, -R4 = —0111 — 30133, 

-R5 = 0222 + 30233, -^6 = "0333 " 30223, Rj = -60123- (34) 

There is a typographical error in eq. (29)3 of Al'shits and Lothe 9j: the term (— 3-R3 — 
2Ri) should be (-3i?3 - R4). 

We mention some other properties of the tensor 0, beginning with 

|0|2 = ir2, where T = (Ai - A2)(A2 - A3)(A3 - Ai) , (35) 



and the related result, from eq. (J2S1), that 

|/(m)|<^|0||m|^ (36) 

with equality iff |a ■ m| = |b ■ m| = |c ■ m| = |m|/\/3. In addition to being traceless, 
eqs. (I23|) and (j24|) imply that (p satisfies the orthogonality-like relations 

0ip<70ipg = T^r^l^ij- (37) 



1^ 



Define the 6- vectors Vj, i = 1, 2, 3, 



T 



Vi = |0ill, 0i22, 0i33, V20i23, V20i31, V20il2| , Z = 1,2,3. (38) 

Then eq. (jTTj) may be rewritten in a form that appears more like an orthogonality 

property: 

I 
V* ■ Vj = -|0|^(5ij, i,j = 1,2,3, (39) 

while the identities ()2fi|l become 

V, ■ r = 0, z = 1, 2, 3, r = (1, 1, 1, 0, 0, 0)^. (40) 

Using for instance, eqs. (fH^ and (jlUj) with i = I, the magnitude of the acoustic axis 
tensor can be expressed by means of no more than five of the ten components of (pijk 
(compare with the expression (plj) that employed 7 coefficients): 

- |0|2 = 02^, + 02^2 + (0111 + 0122)' + 20^23 + 20^,2 + 0^13 • (41) 

The tensor </)(n) must vanish identically if n is an acoustic axes. In principle, this 
requires all seven components c„, n = 0, 1, 2, . . . , 6 must be zero, or all five of the com- 
ponents in eq. ^^- In practice, in order to use the function /(m) to test whether 
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or not two eigenvalues coincide, the direction m must be chosen so that the product 
g{m) = (a ■ m)(b ■ m)(c ■ m) is non-zero. For an arbitrary choice of m this could vanish 
if m lies in one of the planes normal to a, b or c. Given two or three arbitrary directions, 
g could vanish for all if they each lie in these planes. However, for four directions, such 
as ei, e2, e^, and ei + e2 + 63, then at least one of them must yield a non-zero g. Note 
that /(cj) = (piii (no sum), and using identities like 

/(ei + e2 ± 63) = 60123 - 2(0111 + ^222 ± ^333), (42) 

we surmise that the simultaneous vanishing of, e.g., ^m, 0222, ^333 and 0i23, is sufficient 
to ensure the existence of an acoustic axis. 

5 Specific elastic symmetries and weak anisotropy 

5.1 Monoclinic symmetry 

We first provide a general procedure for determining acoustic axes in materials with 
monoclinic symmetry. The off-diagonal elements of Q are, for arbitrary anisotropy, 

Q23 = C^quI + C24.nl + 034^3 + (C23 + 044)^27^3 + (C36 + 045)^17^3 + (C25 + C46)nin2, 
Q3I = Ci^nj + C^quI + Cs^nl + (C45 + 036)^2723 + (Ci3 + 055)721723 + (cm + C56)7Zi7l2(43) 
Q12 = CiquI + C26nl + C45713 + (C25 + C46)7i27l3 + (Cu + €5^)711713 + (Ci2 + C66)7Zl7l2, 

Let 111 be normal to a plane of symmetry, so that the moduli Cj5 and Cj6, i = 1,2,3,4, 
all vanish. It is assumed that the acoustic axis lies in the plane 7i2 = 0, which is possible 
if we allow the original acoustical tensor to be rotated about ei by an as- yet unknown 
angle, say ^o- In this case, 

Q23 = C^gUJ + Cunl, Q31 = (Ci3 + C55)71i7l3, Q12 = (Ci4 + C56)71i7Z3 (44) 

and hence, from eqs. ((Tj), 

72 77/ 

qf = (Ci4 + C56)(C13 + C55) 2^-^ 2' (45) 

C56^i + £34^3 

Q2 = i^^—^)ic56nj + 034713), ql = i^^-—^)ic5enl + 034713), (46) 

Cl3 + C55 Ci4 + C56 

Using Q22 = c^Quf + c^aJiI and Q33 = c^^n\ + 033713 the condition D2 = -D3 becomes 



Ci4 + C56 Ci3 + C55 

C66 ~ C55 + ( ■ 1 ■ ) C56 

Cl3 + C55 Ci4 -+- C56 



71i 



Ci4 + C56 Ci3 + C55 

C33 — C44 + ( ■ ■ ) C34 

Cl3 + C55 Ci4 -\- C56 



%• 



(47) 

Similarly, using Qn = ciin\ + 05572,3, condition Di = D2 and the identity n\ + n\ = 1, 
provides a single but lengthy equation involving the moduli only, which must be solved 
for ^0- 
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5.2 RTHC crystals 

More explicit and useful methods exist for materials with symmetries higher than mon- 
oclinic. Thus, Mozhaev et al. jEl derived a method similar to that for the monoclinic 
material above for trigonal materials, but which gives an explicit expression for sin36'o. 
All other symmetries above monoclinic, viz. orthorhombic, tetragonal, hexagonal and 
cubic, are included in the category of "RHTC crystals". These symmetries have been 
considered in detail by Boulanger and Hayes [H E] who showed, using a result of Fedorov 
and Fedorov |3], that acoustic axes in RTHC crystals can be completely classified by 
relatively simple algebraic means. 

The key to the analysis of Boulanger and Hayes [^ 1^ is the observation of Fedorov 
and Fedorov that the acoustical tensor for RTHC materials can be expressed in the form 
of eq. (0)) with s = +1 where D^, i = 1,2,3 and the vector q have explicit forms: 

3 

A = ^ ttijU^j , q = {PiUi, ^2112, Pans). (48) 

i=i 

The constants aij and Pi can be expressed in terms of the moduli jH |H] . The results of 
Boulanger and Hayes are rederived in Section IHl using geometrical, and perhaps more 
general, arguments which shed new light on the structure of the acoustic axes in RTHC 
crystals. First we consider the separate but important case of materials with weak 
anisotropy. 

5.3 Acoustic axes in weakly anisotropic materials 

A material is defined as weakly anisotropic if the stiffness tensor is of the form C = 
C*^°^ + eC^^\ where C^^^ is isotropic and |e| -C 1. Accordingly, the acoustic axis function 
of eq. (J22I) can be expanded using Q(n) = Q'^°^(n) + eQ*^"'^^(n), as 

[m, Qm, Q^m] = /("^m) + ef^^\ni) + e^f^^\ni) + e^f^^\ni). (49) 

The term f^^^ui) is identically zero by virtue of the fact that C*^'^-* is isotropic, and hence 
the condition for acoustic axes becomes 

/W(m) = -e/(2)(m) - e^/^^Hm). (50) 

After some simplification, it can be shown that 

/«(m) = -(A(°)+/i(°))2(n-m)[m, n, QWPm] 

= (n • m) m • -j/rn , (51) 

where A^*^-* and yU*^"-* are the isotropic Lame moduli, P = I — n C?) n, and the tensor tp is 
defined by 

i^^J = -l{^^'^ + ^^^''^)^{e^klnkS,l + eJMnkSu), and S = PQ^^^F . (52) 
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Equations (fKHjl and (jHT|l imply 

4>ijk = i^ijrik + ipikUj + ipjkrii + 0(e). (53) 

Hence, to leading order in e, the condition for the existence of acoustic axes reduces to 
the requirement that the second order tensor if) vanish. We define the weakly anisotropic 
acoustic axis condition as "j/j = 0. Some insight into tjy follows from the general properties 
of (j), from which it may be shown that i/^n = 0, tr '?/' = and -0 = P tr(i/' )/2. 
Furthermore, it can be shown that i/? vanishes iff S = AP, where A = |trS. Since S is 
rank degenerate by virtue of its definition, ?/? vanishes iff the two non-zero eigenvalues of 
S are equal. Let Ai and A2 be the eigenvalues, then the weakly anisotropic acoustic axis 
condition 2(A^ + A2) — (Ai + A2)^ = can be expressed, using the identities trS" = X^ + X^, 
as 

2trS^ - (trS)2 = 0. (54) 

Equation (j54j) is interesting in that it apparently provides a single condition for 
the existence of acoustic axes. However, its simplicity begs the question of whether 
satisfaction of this weakly anisotropic acoustic axis condition is equivalent to the exact 
condition (j2^- In order to answer this we assume, with no loss in generality, that 
the isotropic medium is chosen such that pc| = n ■ Qn and pc^ = |trPQP, where 
cl = {X^^' + 2fj,^^^Y'^ / p^'"^ and ct = (/i^^Vp)"^ ^^^ the longitudinal and transverse wave 
speeds, respectively. This implies that trS = automatically. The weakly anisotropic 
acoustic axis condition becomes S = 0, and when it is satisfied, it can be shown that 
the acoustic tensor must be of the form Q = Q'^''^ + e(p (g) n + n ® p) for some |p| = 1, 
p ■ n = 0. Thus, pc"^ is an exact eigenvalue with transverse polarization in the direction 
p A n. A simple analysis shows that the quasi-transverse eigenvalue is 

pc'^ = pel - e\pcl - pcl)~^ + O(e^). (55) 

Hence, when the weakly anisotropic acoustic axis condition is satisfied, one of the 
split transverse waves is a pure transverse wave with speed ct, and the other has speed 
Ct + O(e^). The direction n is an acoustic axis iff these speeds are equal, and occurs if 
e = 0, in which case all directions perpendicular to n are transverse directions and a pure 
longitudinal wave exists in direction n. Thus, we conclude that the weakly anisotropic 
acoustic axis condition predicts an acoustic axis iff n is also a longitudinal axis. This 
additional requirement greatly restricts the utility of the weakly anisotropic test for 
acoustic axes, to the point that it is virtually useless in practice. In other words, if 
one calculates slowness surfaces to leading order in the weak anisotropy, the occurrence 
of double roots is not in general an indication of real acoustic axes but coincides with 
acoustic axes only if n is also a longitudinal axis. 
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6 Geometrical interpretation of the acoustic axes 
conditions 

The analysis in Section |2l motivated the previously known form of the conditions for the 
existence of acoustic axes jL . Having found the more general form involving the acoustic 
axis tensor in Sections El and EJ we now return to the ^^ad hoc" analysis based on the 
assumed form (jH)) and show that it also explains the variety of possibilities contained in 
the third order tensor 0. 

Without any loss in generality, we assume that Di > D2 > D^. Then by expressing 
D in biaxial form similar to eq. (^, the acoustic tensor can be written 

Q(n) = D2I + -{d+ d_ + d_ ® d+) + q^ q, (56) 

where 

d± = Dif ei±D^f eg, (57) 

and the abbreviated notation 

D,,=D,-D^, (58) 

is used. Thus, Q is uniaxial when the sum of terms in eq. (jSUIl reduces to the standard 
form for the uniaxial tensor, eq. (j3)). The question can then be phrased as follows: 
assuming Di > D2 > D3, what is the constraint upon q that ensures Q of eq. (j^ is 
uniaxial? The answer requires consideration of various possibilities for q, which we now 
examine in detail. 

First, noting that Di is the largest of the Di, it is clear that there is no way that Q can 
be uniaxial if q || ei. If q is in the 62 (63) direction, then its magnitude must be exactly 
D^2 (D23 ) ill order for Q to be oblate (prolate). However, it is not sufficient that q lie 
along one of the principal directions of D, since it must have a very specific magnitude 

1/2 

if Q is to be uniaxial. For instance, q || 62 and |q| = D^2 implies Qu — Q22 = 
with the addition constraint q || 62, which cannot hold in general, as illustrated by the 
specific example of RTHC crystals, see [5J. The cases of q in coordinate planes and out 
of coordinate planes must therefore be considered separately. If we assume that q lies 
in the plane spanned by ei and 62, then in order that Q is uniaxial we must have 

Diei (g) ei + D2e2 (g) e2 + q ® q = -D' (ei (g) ei + e2 (S) 62), (59) 

for some D'. However, writing 

DiBi ® ei + D2e2 ® 62 

= -{Di + D2){ei ® ei + e2 ® 62) + -{Di - D2){ei (g) ei - e2 ® 62), (60) 

then it it clear from the second term in the right member, specifically the combination 

(ei ® ei — 62 ® 62), that q = giBi + 52^2 cannot satisfy eq. (I^Uj) unless gi = and 
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g| = Di — D2- Thus, we again find that the acoustical tensor becomes uniaxial only if q 
is along a coordinate axes, and then only if a specific relation holds, which cannot hold 
in general. A similar negative result is found if q is sought in the plane spanned by 62 
and 63: the only possible q that leads to uniaxility is q = ±-D23 63. 

We are thus led to consider q in the plane of ei and 63, or equivalently, the plane 
spanned by d_|_ and d_. Assume the form 

for some constants a±, which implies, using eq. fj56|l . that 
1 



Q = ^2l , ^ 



aM+®d+ + ald_ (g)d_ + (d+(g)d_ + d_ (g)d+)(2-a+a„) . (62) 



This reduces to the uniaxial form (^ iff a+a^ = 1, and then Q(n) = D2I + u ® u, 
where u = (a;+d+ + Q;_d_)/2. Thus, at uniaxiality, 

q = 2['^"+~"-)^i2^^i + (a+ + a-)^2f ^3], (63) 

u = l[{a+ + a^)Dli^e^ + {a+-a.)Dli'e,]. (64) 

Removing the dependence upon the parameter «+ (and a_ = l/a+) implies that q and 
u lie on complementary hyperbolae in the 61,63 plane: 



D23 D12 D12 D23 

The eigenvalues are {Dq, D2, -D2}, where the largest eigenvalue Do follows from eq. (jHHjl 
as Do = D2 + u'^ = Di - D2 + Ds + q^. Noting 

g' = (^^ )^Di3 + /^23>^23, " =(^^ )^Di3 + D12 > /^12, (66) 

we see that Do > Di. Also, u = (^3/^12, 0, gi-D23) (-Di2-D23)^^^^, indicating that the 
plane of circular polarization is spanned by 62 and {qiD23, 0, —^3/^12). 

Boulanger and Hayes j3] obtained the condition (jHSjli for RTHC crystals; specifically, 
it is identical to eq. (3.6) of [H]. This is also equivalent to (Qn— Q22)(Q33— Q22)— Q31 = 0, 
which in turn is equivalent to 0ii2 = or ^332 = 0. Equation (jHKjl i is an instance of the 
specific application of the condition (puj (no sum), as compared with the original derived 
requirements (lllj). 

We now consider q out of the coordinate planes, so that no component is zero. The 
eigenvalues of Q are the zeroes of g{X) = det(D + q q — AI). This cubic in A takes 
on the values g{Di) = qfDijDik, where {i,j,k} is any permutation of {1,2,3}. Thus, 
g{Di) > 0, g{D2) < and g{D^) > 0, indicating that the eigenvalues of Q must satisfy 
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D3 < A3 < 1^2 < -^2 < -Di < Ai. In particular, the eigenvalues are distinct as long as 
?i?2?3 7^ and Di > D2 > D3 (see Section 4 of 4]). Therefore, the latter constraint 
must be relaxed and either D12 = or D23 = assumed, that is 

Q = D2I + ^12 ei ® ei + q O q or Q = D2I - /^23 63 ® 63 + q O q, (67) 

respectively. The only way (jEZJi can be uniaxial is if gi = and g^ = D12, implying the 
oblate Q = Dil — D^e ® e where e = ei A q/|q|. Similarly, in order for (jUTj) ^ to be 
uniaxial, q must be aligned with e^. Both solutions again place q in coordinate planes, 
and therefore the only possibility for uniaxiality with q out of the coordinate planes is 
that D is spherical, that is, Di = D2 = -D3. This was derived previously in Section |21 
and was also found by Boulanger and Hayes for RTHC crystals, see eq. (4.7) of 5j. 

Finally, we consider the possibility q = 0, in which case Q is diagonal and uniaxiality 
occurs when two of the diagonal elements are equal. As noted in Section |5 the only 
nonzero element of is 0i23 when Q is diagonal. Hence, this particular case of q for 
acoustic axes corresponds to the unique condition (pus = 0. 

6.1 Discussion 

We have seen that the geometrical derivation of the acoustic axes conditions reproduces 
all elements of the acoustic axis tensor 0. These results provide new insights into the 
structure of symmetric second order tensors, beyond the specific interest in degenerate 
acoustical tensors. Consider, for instance, the stress tensor of a body in a state of 
triaxial stress, for which the principal stresses are distinct, o"i > (72 > 0-3. The body 
is then subject to an additional uniaxial stress of the form Tn ® n, which is tensile 
or compressive depending as T is positive or negative, respectively. Suppose that the 
values of ai, (T2, (J3 and T are all fixed (given), and the question is posed: what choice or 
choices of the direction of uniaxial additional stress will result in a state of total stress 
that is uniaxial? The answer follows immediately from the acoustic axes analysis: the 
direction n must lie in the plane of the least and greatest stress, i.e. n = riiei + n^e^, 
where, by analogy with eq. (jU^ i. 

"^ "^ ' (68) 



0"2 — cr-i cji — a2 T 

The state of uniaxial stress can only be achieved if T > a2 — as or T < — (cti — 0-2), 
and the resultant stress is 

cr = cr2l + (ctq - 0-2) e (g) e , (69) 

where the third stress ctq and its direction e are given by 

rj.. , „ (ai - aa) ^1361 + (^2 - ^3) ^163 .„„^ 

CTq = i + (Ti - 0-2 + CTS, e = sn 2 I / \2 211/2 " ^'^) 

[{ai - 0-2)^^1 + ((T2 - asynfY/^ 
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It follows that the stress is prolate (oblate) with ctq > o"! {o'o < os), if the applied 
uniaxial stress is tensile (compressive). Note that the stress analog problem is slightly 
different from the acoustic axes issue since T can be negative as well as positive. 

7 Conclusion 

Conditions for the existence of acoustic axes can be formulated in terms of a simple 
geometrical framework, as in Sections |21 and (jH)), or in terms of the rigorous requirements 
on the third order acoustic axis tensor (f) of SectionHJ Specifically, </> = at acoustic axes, 
which in turn involves a multiplicity of identities some of which are often automatically 
satisfied, as in symmetry planes. The case of RTHC crystals is best understood in terms 
of the geometrical approach. 

Appendix 

A Proof of equation (^ 

The unit vectors ai and a2 are determined by finding isotropic bivectors vi, V2 that solve 

V ■ Av = 0. (71) 

A bivector [8J is a complex valued vector, v = p + iq. The dot product is defined as for 
real vectors (no complex conjugates are involved), thus, Vi ■ V2 = Pi ■ P2 — Qi ■ q2 + *(Pi • 
q2 + P2 ■ Qi)- The bivector v = p + iq is called isotropic if v ■ v = 0, and occurs iff p 
and q are orthogonal and of equal magnitude. Zeroes of eq. (f7T|) may be sought in the 
form V = (z^ — l)ei + i(z^ + l)e2 + 2ze^, which is clearly an isotropic bivector for all 
complex-valued z. Equation ()7H1 is a quadratic in v, homogeneous in the components 
of V, and becomes a quartic in the complex variable z. The roots occur in pairs, since it 
is clear that if v is a zero, and so also is v*, hence —1/z*, which gives the same bivector. 
Thus, two distinct bivectors, {vi, V2}, solve eq. (f7T|) . although they are not unique 
since any other pair, {v'^^, V2} = {e^'^^vi, €^'^'^^2}, ipi, 'ip2 £ R, also solves eq. (f7T|) . The 
non-uniqueness indicates that the planes spanned by the bivectors are of importance, 
rather than the bivectors per se. The real unit vectors ai, a2 are the normals to these 
planes, since this choice of directions for ai, a2 ensures that A of eq. (^ has vi, V2 as 
zeroes in eq. (f7T|) . It may be shown that cosh7j a^ = cos6'j ei -|- sin6'j 62 + sinh7j 63, 
where Zj = e~*^'^J"'"*^^\ j = 1, 2. Let d = ai A a2/|ai A a2|, and referring to eq. (H)), 

A = d-Ad, A' = (tr A- 3d- Ad)/ai-a2 . (72) 

This completes the construction of the general form of A in eq. (^. 
It is interesting to note that the bivectors vi and V2 satisify 

Avj = Avj + fij^j no sum, (73) 
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where /xi = (A'/2)a2-vi and /i2 = (A'/2)ai-V2. With no loss in generahty, we may express 
the bivectors as Vj = e*'^^(d + zd Aa^) (no sum) for some ipi, ip2 & R- A straightforward 
calculation gives 

^^ = ie'^'{X'/2)\sii A aal, fi2 = -le'^f^^iX' /2)\sii A as], (74) 

indicating that |/ii| = |yU2|- If V'l = "^2 ± tt = ■?/' then /Ui = /i2 = yU and we have the 
general result 

Avj = \\j + ^Rj where /i = ie'''^{X' /2)\ai A a2|. (75) 

In particular, the bivectors are eigenvectors with the same eigenvalue iff the bivectors 
are coplanar, i.e. iff ai A a2 = 0. This reproduces the result of Boulanger and Hayes |H] 
that A is uniaxial iff the the bivectors are also eigenvectors. 

B Cartan decomposition of the acoustic axis tensor 

A procedure is described for determining the Cartan decomposition of the tensor 0. 
The method, which can be generalized to tensors of other order, is based on finding the 
irreducible representation of the elements of (f) under the action of SO (2). Let 63 be 
the preferred direction for rotation under S0(2), and write m = maGa + fn^e^, where 
a = 1,2. Thus, 

/(m) = (pap-y m^mpm^ + 30c,/33 m^mpm^ + 30^33 m^ml + ^333 m^, (76) 

where the fact that (f) is totally symmetric has been used. Introduce the angle 9 via 
mi = r cos 6, m2 = r sin 6, where r^ = 1 — ml, then 

0a/37 m^ampm,^ = (0111 cos^ 9 + 0222 sin^ + 30112 cos^ 9 sin 9 + 30i22 cos 9 sin^ 9) r^, 
30a/33 "^a'^/3'^3 = (0113 cos^ + 0223 sin^^ + 20123 COS ^ sin 6')3r^m3 (77) 

30a33 m^ml = (0133 COS 9 + 0233 siu 9) 3rml . 

Applying simple trigonometric relations, implies 

/(m) = -(0111 -30122)?^^ cos 36^ + -(-0222 + 30ii2)r^ sin 36^ 

3 

+ (0113 - 4>22z)^r^m^ cos 29 + 0i233r^m3 sin 29 

+ ["^30133 + Y^*^!" + 0122)]3rCOs6' + [m30233 + ^(^^222 + 0112)]3r slu 6' 

3 

+ (0113 + 0223)2^^"^3 + 0333"^i- (78) 

Using the harmonic property (|^. it can be deduced by inspection of eq. (fTHj) that the 
quantities 

0333) 
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(f)l33 ± i(j)233, 

0113 -0223±»i2(/)i23, 

0111 - 30122 =F «(0222 - 30112), (79) 

correspond to the coefficients of the partition of the 7— dimensional space of third order 
totally symmetric harmonic tensors into 1 + 2 + 2 + 2 subspaces as in eq. (P7j) . 
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